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Cosmic Neutrino Last Scattering Surface 
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Neutrinos decoupled from the rest of the cosmic plasma when the Universe was less than one 
second old, far earlier than the photons, which decoupled at t = 380, 000 years. Surprisingly, though, 
the last scattering surface of massive neutrinos is much closer to us than that of the photons. Here 
we calculate the properties of the last scattering surfaces of the three species of neutrinos. 

PACS numbers: 95.35.+d; 95.85.Pw 



I. INTRODUCTION 

The standard cosmological model predicts that neutri- 
nos were produced in the early universe and are present 
today with an abundance 1 of 112 cm -3 per species 
Detecting this background remains a tantalizing experi- 
mental dream d, [1,0, d, H> fill- w ith recent develop- 
ments encouraging the optimists [l2l [l3l . [bH . [TBI . [l6l . u7\ , 

m. 

For the time being, the question of where the cosmic 
neutrinos come from remains an academic question. Yet 
it is of sufficient interest that, even if there were no chance 
of detection, the origin of the cosmic neutrino background 
(CNB) seems worthy of theoretical study. 

The neutrino Last Scattering Surface (LSS) is typically 
thought of as being located a given distance from us with 
a small but finite width, similar to the last scattering sur- 
face of the cosmic microwave background (CMB). Neu- 
trinos last scatter when the temperature of the universe 
was a few MeV and the universe was less than a second 
old, while the photons in the CMB last scattered much 
later when the temperature was 1/3 eV at t = 380,000 
years, so it is natural to assume that the neutrino LSS 
is further away than that of the CMB. Calculating how 
much further away leads to a little surprise. For, even 
massless particles can travel only a very small comoving 
distance in the very early universe, so the CMB comes 
from a comoving distance about 9540 h~ x Mpc away from 
us (where the present expansion rate is H = 100 h km 
s -1 Mpc -1 ), while massless neutrinos arrive from a co- 
moving distance 9735 h~ L Mpc away. That is, neutrinos 
travel only about 200 Mpc (comoving) in the first 380,000 
years. 

We show in this paper that this slightly surprising re- 
sult evolves into another counter-intuitive result for mas- 



1 The uncertainty on this prediction is set by the sub-percent un- 
certainty in the temperature of the cosmic microwave background 
temperature pj, which serves to calibrate the thermal number 
density. The predicted abundance is roughly one percent larger 
than thermal due to e+/e _ heating 



sive neutrinos [19|, [20( : the CNB actually reaches us from 
closer than the CMB. Even for neutrino masses as small 
as 0.05 eV (and one of the neutrinos must be at least this 
massive) the effect is dramatic. For neutrinos with mass 
of 1 eV, the effect is truly striking with most of these 
neutrinos arriving from only several hundred Mpc away! 



II. LAST SCATTERING SURFACE OF 
MASSIVE NEUTRINOS 

Neutrinos stopped scattering when temperatures were 
of order a few MeV and the universe was less than a 
second old. At that time each neutrino species had a 
Fermi-Dirac distribution for a massless particle (assum- 
ing - as we will throughout - zero chemical potential and 
no heating from electron-positron annihilation). 

Consider first the calculation of the last scattering sur- 
face of a massless neutrino. The comoving distance trav- 
elled by a massless particle starting from i, = 1 sec until 
today is 
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Note that early on H(a) = HqQ^^ a" 2 where il r = 
8.3 x 10~ 5 is the radiation density today in units of the 
critical density. The integrand peaks at late times, so the 
contribution to the comoving distance from early times 
(fii ~ 1 sec) is negligible; this explains our nonchalance 
in defining the initial time. It also explains why all three 
species of neutrinos share the same last scattering surface 
even though electron neutrinos decouple slightly later 
than do the other two species. The comoving distance 
travelled by neutrinos until the time of photon last scat- 



tering at a, 



1090.5 ±0.95 \2M is: 



1 



X* = 
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y/a + a cq 



(2) 



where VL m is the matter density today in units of the crit- 
ical density and a eq is the value of the scale factor when 
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Using 
this equates 



the densities of matter and radiation are equal 
the standard cosmological parameters y 
to the 200 hr 1 Mpc difference alluded to in the introduc- 
tion. 

Massive neutrinos slow down once they become non- 
relativistic, so the integral determining the distance to 
the last scattering surface generalizes to: 



junction: 



dt 



Pol a 



u ">(*) V(Po/ a ) 2 + m l 



(3) 



where the second term in the integrand is the redshifted 
velocity p/E, with po the current neutrino momentum. 
The neutrino temperature today is T v — 1.95 x 10~ 4 
eV, so there will be a range of pn's drawn from a Fermi- 
Dirac distribution, each of which will be associated with 
a different distance to the LSS. Fig.[T]plots this distance 
as a function of neutrino mass for two different values of 
the present day neutrino momentum. 
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FIG. 1: The comoving distance travelled by a massive neu- 
trino since decoupling until today as a function of mass for 
two different values of the neutrino momentum. Neutrinos 
with smaller momenta are travelling more slowly and there- 
fore travel a shorter distance. Note that for masses above 
10 -4 eV, the neutrino LSS is much closer than that of the 
CMB. 

Since neutrinos with different momenta arrive from dif- 
ferent distances 2 , the last scattering surface of the CNB 
is quite broad compared to that of the CMB. To quan- 
tify this, we can define the probability that a neutrino 
last scattered a distance \ away from us, or the visibility 



2 Ref. [22l | also mentioned this feature of neutrinos and proposed 
to exploit it to test the Copernican Principle. 



dP _ dP ( d\ 
dx dp \dp 



(4) 



where the equality uses the chain rule; the first differ- 
ential probability is given by the massless Fermi-Dirac 
distribution 12311 : 



Pi 



dP 



d Po 3C(3)T8ePo/ T " + l' 



(5) 



and the second term on the right is obtained by differen- 
tiating Eq. ©. 

Fig. [2] displays the probability that a neutrino with 
a given mass arrived from a distance Note for all 
masses above 10 -4 eV, the spread in arrival distances is 
much larger than the spread in the CMB last scattering 
surface. The width of the CNB and CMB last scattering 
surfaces have different origins: the CMB does not have an 
infinitely thin last scattering surface because the process 
of recombination, and therefore decoupling, extends over 
a finite time period. The CNB last scattering surface 
is thick, reflecting the different velocities of the different 
momenta in the Fermi-Dirac distributions. 
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FIG. 2: The probability that a neutrino with mass m last 
scatters at a given comoving distance from us (the visibility 
function). Massive neutrinos travel more slowly than massless 
neutrinos so arrive here from much closer distances. Also 
shown is the last scattering surface of the cosmic microwave 
background, virtually indistinguishable from that of an m v = 
10~ 4 eV neutrino. 

One might ask about the last scattering surface of 
heavier particles, such as sterile neutrinos or ordinary 
cold dark matter particles. Those that become non- 
relativistic before equality travel a distance of order the 
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comoving horizon at equality (136 Mpc) times the co- 
moving velocity at equality [{p / m) / a-^q\. A keV sterile 
neutrino therefore would have a last scattering surface of 
order a Mpc away. For particles more massive than this 
- e.g., Cold Dark Matter - their LSS is so close that the 
distance traveled in a straight line (before gravity in our 
halo started moving them around) was negligible. So in 
some sense, the question ceases to make much sense for 
masses above a keV. Some of these considerations apply 
even to light neutrinos in our Galaxy: to determine what 
has happened to these neutrinos recently, one must c arr 
out simulations along the lines of those presented in 



III. OSCILLATIONS 



arry 

m 



Neutrinos are produced in flavor eigenstates but then 
propagate as mass eigenstates. Therefore, the simple no- 
tion of a neutrino with fixed mass having its own last 
scattering surface is a bit too naive. As emphasized in 
[25l | . the number density of neutrinos in mass eigenstate 
i is: 



drij = 



a=l 



\U ai \ dn D 



(6) 



where U a i are the elements of the unitary matrix which 
transforms from the mass basis to the flavor basis, and 
the sum is over all three flavors. 

In principle the number densities of the three neutri- 
nos flavors dn a could differ if, e.g., each had a different 
chemical potential [2a]. Here we assume that the chem- 
ical potentials are very small so all the dn a are for all 
practical purposes identical. In that case, they can be 
removed from the sum, and then the sum over the uni- 
tary matrix elements squared simply gives unity. So 



drij = dn n 



dp 



2tt 2 ePo/T„ + 1 



(7) 



independent of mass eigenstate. So we expect the calcu- 
lation of the previous section to reflect the different last 
scattering surfaces of the different mass eigenstates. 

Detection, however, will take place in flavor space, with 
electron neutrinos and anti-neutrinos. To compute the 
last scattering surface of these eigenstates, we must trans- 
form back to flavor space, weighting each mass eigenstate 
by its own visibility function: 



dP, 

dx 



dP 
dpo 
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Fig. [3] shows the resulting probability for an electron 
neutrino with two possible mass schemes assuming the 
tri-bimaximal mixing matrix [27], [28| (U e ± — 
U e 2 — 1/v3j and U e3 — 0). Note the interesting dou- 
ble peaked structure in the normal case, a signature of 
the quantum mechanical oscillations which dictate that 



there is roughly a 2/3 probability the electron neutrino 
propagates with very small mass mi and 1/3 with mass 
m 2 . 
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FIG. 3: Probability for the distance to the last scattering 
surface for an electron neutrino with normal hierarchy (mi = 
0.009 eV; mi = 10 -4 eV; sin #13 = 0) and degenerate masses 
(mi = mi = m.3 = 0.2 eV). 



IV. CONCLUSIONS 

The last scattering surface of the cosmic neutrino back- 
ground is much broader and much closer than that of the 
cosmic microwave background. Indeed, as depicted in 
Fig. [3l the last scattering "surface" of the electron neu- 
trinos and anti-neutrinos has a very rich structure due to 
oscillations. 

Are there any observable consequences of these dis- 
tances to the last scattering surface? We can think of 
three potentially interesting follow-ups: 

• Neutrino Acoustic Oscillations 

The well-known Baryon Acoustic Oscillations 
(BAO) arise because the photon-baryon gas trav- 
els a distance equal to the sound horizon at de- 
coupling, after which the baryons stop. In real 
space, this leads to an overdensity of baryons in 
a spherical shell surrounding an initial overdensity, 
and ultimately to a bump in the correlation func- 
tion. Massive neutrinos would seem to share many 
of the same features: adiabatic perturbations lead 
to neutrinos initially being overdense at the same 
places as the baryons; neutrinos travel a finite dis- 
tance since last scattering; and this distance might 
show up as a feature in the power spectrum. Unfor- 
tunately, we see no evidence for this feature when 
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running the Boltzmann code [2!| which solves for 
the linear evolution of perturbations (and by track- 
ing different neutrino momenta implicitly encodes 
all of the relevant physics). The difference between 
neutrinos and the baryons are that the neutrino 
last scattering surface is very broad due to the 
Fermi-Dirac distribution; this breadth smooths out 
any (very small) feature rendering it undetectable. 
Searching for a feature might still be useful as a way 
of constraining the distribution function of neutri- 
nos in the CNB. 

• Galactic Distribution function 

For direct detection, it is important to know the 
distribution function of neutrinos in a Milky- Way 
sized galaxy. Are neutrinos overdense? Has their 
momentum distribution changed due to virializa- 
tion? If capture has occurred, the distance to 
the LSS would likewise be modified. Indeed, it is 
possible that there would be a directional effect: 
neutrinos from the center of the Galaxy might be 
captured and hence have relatively small distances 
to the LSS, while those from the poles maintain- 
ing their primordial distribution functions and dis- 
tances to LSS. 

• Anisotropics 

If experimentalists do succeed in detecting the 



CNB, in the far future one might imagine maps 
of the anisotropy of this background [30(. If the 
neutrinos have been sloshing in our Galaxy for a 
number of orbits, this angular information might be 
lost. But some of the neutrinos - the lightest ones 
or those arriving away from the Galactic Center - 
will arrive undeflected (see, e.g., the right panel in 
Fig. 2 in Ref. [31J which suggests very little trap- 
ping of neutrinos with masses less than 1 eV in a 
Galaxy of our size). These would provide informa- 
tion about overdensities at a very early time from 
locations about which we will already have much 
information from galaxy surveys. The galaxy sur- 
veys probe the same locations but at much later 
times. This would be an almost unique opportu- 
nity to probe the evolution of structure in given 
regions. 
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